An Order Parameter for the Mott-Hubbard Transition in One Dimension 
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We propose an order parameter for the one dimensional Mott-Hubbard transition and provide nu- 
merical evidence and general theoretical arguments for the correctness of our proposal. In addition, 
we discuss some of the implications of this picture of the one dimensional Mott-Hubbard transition 
and speculate about possible higher dimensional analogs. 
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The one dimensional Hubbard model: 

Hid Hubb = t ^2 c i,x+v c cr,x + n h^ n i,x (1) 

x.a.u— ±1 cc 

is known to exhibit a T = quantum transition be- 
tween metal and insulator as a repulsive interaction in 
the model is turned on || and an analogous transition is 
believed to occur in higher dimensions |2j. These "tran- 
sitions" are somewhat peculiar in that a physical, local 
order parameter has not been identified. Instead, fol- 
lowing Kohn's suggestion ||, the "order parameter" for 
these transitions is generally taken to be the weight in 
the zero frequency delta function in the electrical con- 
ductivity at zero temperature. In the metallic phase this 
weight is finite and in the insulator it vanishes, so that 
the proposed order parameter does in fact distinguish 
the two phases qualitatively as it should. However, the 
conductivity is clearly not locally defined, so that it is 
meaningless to speak of, for example, a diverging length 
scale in its correlation functions, and it leads to a qual- 
itative distinction between metal and insulator only at 
zero temperature. In the limit of large spatial dimen- 
sions, the Hubbard model undergoes a finite tempera- 
ture metal-insulator transition associated with a crossing 
of the free energies of two distinct phases 0]; however, 
there is no order parameter distinguishing the phases in 
that case, and the transition line vanishes at a critical 
point, much like the liquid-gas transition of water. This 
is not expected to be the case for the metal insulator 
transition in one spatial dimension and may also not be 
generic for finite dimensions. Therefore, it would be de- 
sirable to define some local operator whose expectation 
value or correlations clearly distinguished the metallic 
and insulating phases; the purpose of this paper is to 
propose an operator which fulfills these requirements for 
the one-dimensional transition and to examine possible 
implications for higher dimensions. 

The operator we propose is an extension of the opera- 
tor introduced by us and P. W. Anderson || as exhibiting 



off-diagonal long range order (ODLRO) in a wide vari- 
ety of quantum spin chains. There we noted that the 
groundstate of a generic, gapless XX Z model on N + 2 
sites has a finite overlap in the thermodynamic limit with 
the state obtained by taking the groundstate of the same 
model on N sites and tacking a singlet pair of spins onto 
the end. If we define an operator, 0'(i), which adds a 
pair of sites with their spins in a singlet configuration at 
site i into spin model, and its conjugate operator, 0{i), 
which removes a pair of sites, i and from the model if 
they are in a singlet configuration or else annihilates the 
wavefunction on which it acts, then, as a consequence of 
the overlap mentioned above, this operator has ODLRO 



in its correlation functions: lim 



,<Ot(i)0(j)> * o 



However, limu_ J -i_ >00 (0'!'(i)0( i 7)) ^ C\\ rather, since 
the groundstates of generic XX Z models on N sites and 
N + 2 sites have momenta differing by ir, (O^ (i)O(j)} ~ 
(— 1) % ~ 3 and our order parameter exhibits a broken Z2 
symmetry. This Z2 is unconnected with the division of 
a Neel ordered state on a bipartite lattice into sublat- 
tices: there is no Neel order in one dimension even at 
zero temperature Q. However, the theorems of Mermin 
and Wagner and Coleman [Q do allow the breaking of a 
discrete symmetry, such as the Zi breaking indicated by 
our order parameter. 

What, then, is the nature of this broken symmetry? 
In [9, we utilized the connection between the Heisen- 
berg model and the large U Hubbard model to moti- 
vate the conjecture that the insertion of a single spin 
of spin species a into a Heisenberg model at site j was 
equivalent in the Luttinger liquid description of that 
models low energy physics to the action of the operator 
i j exp[i6 fl!(T (ja)] + (-i)J 1 exp[|e i)(T (ja)]. This not only 
accounts for the ODLRO of the singlet insertion, but ex- 
plains the behavior of the family of operators in which 
singlet pairs of spins are inserted into sites which are 
not adjacent. We have further investigated this proposal 
using exact diagonalization and Haldane H and Shas- 
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try's H solution of the inverse squared exchange model, 
to show that spin insertion is essentially equivalent to 
spinon creation |ic| or the annihilation of a spinon of 
the opposite spin species and the conjecture is now firmly 
supported. However, it leads to a rather peculiar picture 
for the cause of the broken symmetry. The "operator" 
responsible in Luttinger liquid description is (— 1) J , i.e. a 
c number [jll| . This suggests that the order is not the or- 
der of a spin model at all, and we propose that it should 
rather be thought of as the charge "order" which sets in 
when a Mott-Hubbard transition occurs in one dimension 
and the low energy effective theory is transformed from 
one of interacting fermions to one containing only spin 
degrees of freedom. In support of this, we note that all 
of the spin models for which we have demonstrated this 
order can be obtained as the low energy limits of fermion 
models which have undergone such a transition. 

To check the proposal that the order in question is 
really in the charge sector, we have carried out several 
numerical and analytic investigations. First, we have ver- 
ified that the "order" probed by the singlet insertion is 
not that of the spin degrees of freedom by conducting 
a finite temperature study of the singlet insertion-singlet 
deletion correlation function for the XY model. In that 
case, the mapping of the model onto free fermions allows 
us to study large systems sizes at finite temperatures to 
examine how the spin and Zi orders are degraded. Some 
results from Monte Carlo are shown in Fig. |l|. It is clear 
that the order is not degraded by finite temperature in 
the same way as the more conventional spin correlations 
are. 
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FIG. 1. Zero and finite temperature results for the singlet 
insertion/deletion correlator and the more conventional two 
point function of s x . Results are for a one dimensional XY 
model of 100 sites with periodic boundary conditions (PBC). 
Inset: results for correlators for a fixed separation of 25 sites 
as a function of temperature, T. 

Apparently, the identification of the singlet insertion 
with a c number in the spin model does not hold rigor- 
ously and the singlet-correlation/singlet-deletion correla- 
tor is affected by the finite temperature in the spin model, 



lacking ODLRO at high and probably all finite temper- 
ature. Nonetheless, the degradation is enormously less 
severe than that occurring for other, typical spin correla- 
tion functions. For example, at T = .15J, (s x s x ) decays 
exponentially with a correlation length of ~ 4.0 lattice 
spacings, while (O t O) decays by well less than a factor 
of 2 over a separation of 50 sites! It is not even clear 
from results at this temperature that the ODLRO has 
been destroyed. We attribute this difference in behavior 
from conventional spin correlators as the result of the fact 
that the XY model, viewed as a model with an Mott- 
Hubbard ordered charge sector, has an infinitely large 
charge gap. Thus the charge order, as probed by our sin- 
glet insertion, is unaffected by finite temperature, which 
rapidly destroys all of the usual spin order; the singlet 
insertion/deletion correlator is only indirectly connected 
with the charge sector order and it appears that only if 
the spin degrees of freedom are effectively at zero tem- 
perature does it exhibit ODLRO, so its correlator docs 
decay at finite temperatures, but extremely slowly. 

Having shown that our order is almost not degraded 
for separations and temperatures where spin order is to- 
tally lacking, it remains to show that it is sensitive to 
charge ordering of the Mott-Hubbard type. First, it is 
possible to show that, for free, SU(2), spin \ fermions at 
half filling, the operator which adds two adjacent, occu- 
pied sites with spins in a singlet configuration does not 
have ODLRO. Rather, the correlation function decays al- 
gebraically like x~ 1 (or, for SU(N) generalizations, like 
for N inserted sites in an SU(N) singlet). We have 
shown this both numerically and analytically and 
therefore know that free fermion systems, which should 
not show Mott-Hubbard order, do not show ODLRO as 
probed by the singlet insertion operator. 
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FIG. 2. Singlet insertion/deletion correlators for SU(2) 
and SU(3) Hubbard models of various sizes (PBC). Insertion 
and deletion were separated by half the system size, L. 

Conversely, if we consider the large U limit of the one 
dimensional Hubbard model at half filling, then the sin- 
glet insertion operator clearly does display long range 
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order. This can be seen from the Ogata-Shiba factor- 
ized wavefunction |l3|] which has a charge wavefunction 
given by that of spinless fermions and a spin part given by 
the groundstate of the Heisenberg model with the "sites" 
defined as the positions of the spinless fermions. At half 
filling, there is a fermion on every site and the charge part 
is trivial, a property unchanged by inserting or deleting 
singly occupied sites. For this reason the singlet inser- 
tion correlations, are identically those of the Heisenberg 
model and have ODLRO, as we have previously estab- 
lished ||. For finite U, the Lieb-Wu solution Q shows 
there is a charge gap for arbitrarily small U and thus 
one expects Mott-Hubbard order for any positive U, but 
with a magnitude vanishing as U — > 0. This is consistent 
with our numerical results (Fig. ||) for the behavior of 
the singlet insertion/deletion correlator. 

In contrast, for the simplest SU(3) generalization of 
the Hubbard model, a finite critical interaction strength, 
U c , should be required to induce a metal- insulator tran- 
sition [Q . Our numerical results for this model are con- 
sistent with the absence of ODLRO for SU(3) singlet 
insertion for U < U c ~ 2t, and the presence of ODLRO 
for SU (3) singlet insertion for larger values of U. 

We can also make contact between our singlet insertion 
operator and the bosonization picture of the one dimen- 
sional Mott-Hubbard transition; in the large U limit of 
the SU(2) Hubbard model, the singlet insertion opera- 
tor at site j is effectively acting like a constant times 

(— 1)^*<j ' " fc , where is the number of fermions at 
site k. In the Luttinger liquid approach to the Hub- 
bard model, (— l)^-' fc <5 nk maps to cos(9 Pi j)(ja). The 
two point, equal time correlation function of this opera- 
tor decays asymptotically like x" 1 for free particles, as 
does the singlet-insertion singlet-deletion correlator, and, 
in the Luttinger liquid approach to the Hubbard model, 
cos(Op j) takes on an expectation value oc (— 1) J for any 
finite U, with the magnitude vanishing as U — ► via 
a Kosterlitz-Thouless transition. This further supports 
the singlet insertion operator acting as an order parame- 
ter for the one dimensional Mott-Hubbard transition and 
we feel confident that our operator is an effective probe 
of the Mott-Hubbard transition in one dimension. 

What does this picture imply for the Mott-Hubbard 
transition in higher dimensions? First, note that the sin- 
glet insertion operator reveals the presence of a broken 
Z2, spatial symmetry. Such a broken symmetry dou- 
bles the unit cell and might be expected to lead, for 
a half-filled system, to a filled band and a gap to all 
excitations; however, if the broken symmetry involves 
only the charge degrees of freedom, then the result is 
an anti-ferromagnetic insulator rather than a band in- 
sulator, which would have gapped spin degrees of free- 
dom. This behavior is reminiscent of the formation a 
spin density wave (SDW) insulator, however, in that case 
the "order parameter" transfers real electrons across the 



Fermi surface and therefore couples to both spin and 
charge degrees of freedom. Further, the low lying spin 
degrees of freedom of the SDW are Goldstone bosons re- 
sulting from the broken SU(2) symmetry and therefore 
the spin-density wave picture can not apply in one spa- 
tial dimension or in two dimensions at any finite temper- 
ature Q . The breaking of a discrete symmetry, such as 
Z2, is allowed in one dimension at T — and causes the 
Mott-Hubbard transition. Such a breaking is also allowed 
at finite temperature in two or more dimensions. This 
suggests that higher dimensional Mott-Hubbard transi- 
tions in pure Hubbard models may well be finite tem- 
perature transitions at which a Z2 symmetry breaks in 
the charge sector, leaving a low energy effective theory 
containing only magnetic degrees of freedom. This con- 
trasts sharply with the conventional wisdom on the Mott- 
Hubbard transition which follows the line of argument 
that, if the only order parameter for the Mott-Hubbard 
transition is the conductance, then no true, finite tem- 
perature distinction exists between metal and insulator 
|fl5f . This argument fails if one is prepared to admit the 
possibility of the breaking of a Zi symmetry only in the 
charge sector. Based on our one dimensional result that 
the Mott-Hubbard transition does admit a pure charge 
order parameter which breaks a discrete symmetry, we 
feel that this proposal should be seriously considered. 

Interestingly, two organic conductors, (TMTTF) 2 SbF 6 
and (TMTTF) 2 Re0 4 @ exhibit abrupt changes in 
their charge properties consistent with finite temperature 
metal-insulator transitions, and, moreover, these tran- 
sitions are unaccompanied by any detectable magnetic, 
structural or charge density wave transitions. Existing 
explanations for the transition between metallic and in- 
sulating behavior in similar materials |l7j are based on a 
finite temperature crossover that can not account for such 
sharp "transitions" , thus it appears very natural to in- 
terpret this behavior as finite temperature Mott-Hubbard 
transitions of the type we are proposing. 

To test whether the exotic sort of ordering we are 
proposing actually occurs requires a suitable probe for 
Mott-Hubbard order beyond one dimension. In this case, 
the singlet insertion operator we have used is no longer 
appropriate, however, it can be generalized into a useful 
probe. Consider the two dimensional case of a Hubbard 
model of size L x x L y : 

H2D Hubb — y ] ^\\ c lr,x+u,y c o-.x,y + t ^o.x ,y+v c a ,x ,y (2) 
x,y ,(T,y 

+ u z2 

x,y 

where we identify x = L x + 1 with x = 1, with L x an even 
integer, and y = L y + 1 with y = 1. Now, we define a 
multi-singlet insertion operator which inserts L y singlet 
pairs, all at fixed x, and study the dependence on L x of 
the expectation value of this operator acting at x = and 
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the conjugate removal operator acting at x — , 
system has completely anisotropic hopping (tj_ 



L * If the 
0) so 



that each of the L y chains is independent then, at zero 
temperature and for positive U , this expectation value 
behaves like the L y th power of the one dimensional sin- 
glet insertion correlator. For large L x and odd L y , this 
depends on L x only as (— l) - ^, reflecting the underly- 
ing one-dimensional Mott-Hubbard order. Now imagine 
that we turn on some finite interchain hopping, t±, with 
U S> in and U 3> t±. We can use the multi-singlet inser- 
tion to see if the Mott-Hubbard order persists at finite 
t±. 

As a first step in this direction, we considered multi- 
singlet insertion for coupled spin chains; the low energy 
theory for the two dimensional Hubbard model at suffi- 
ciently large U should be that of the Heisenberg model: 



(3) 



with J|| ~ |l and Jj_ ~ |^ 



_l ~ -ffi. The interchain spin cou- 
plings are relevant operators |l8| and could destroy the 
singlet insertion order; however, recall that our order pa- 
rameter acted as a c-numbcr in the spin sector. In this 
case, the usual arguments about perturbed conformal 
field theories imply that its behavior should show no in- 
frared divergences caused by the relevant interchain cou- 
plings (p]]. Therefore the behavior of our multi-singlet 
insertion operator should be be qualitatively the same as 
the isolated chains case: oc (— 1) 2 . If so, the two di- 
mensional Heisenberg model, and, by inference, presum- 
ably the two dimensional Hubbard model in the Mott 
insulating phase, possess hidden, broken Z2 symmetries. 
We have numerically studied the behavior of the multi- 
singlet insertion correlator for two and three leg spin 
chains with the results depicted in Fig. [| 
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FIG. 3. Results of multi-singlet insertion/deletion cor- 
relator for two and three leg spin ladders (PBC, and 
Jout of chain = 2J in chain). Three leg ladder shows alternation, 
and both two and three leg ladders appear to have ODLRO 
for O with the magnitude of (O^O) only slightly reduced from 
that of uncoupled chains. Contrast this with the dramatic ef- 
fects of the interchain coupling on (s z s 2 ), which, for example, 
switches from algebraic to exponential decay for the two leg 
ladder. 

As can be seen from the figure, the relevance of the 
interchain coupling does not lead to any strong change 
in behavior in the multi-singlct insertion operator. For 
the three leg case a broken Z 2 appears to persist; this 
suggests that such a discrete symmetry breaking may be 
present in the two dimensional Hubbard model and un- 
derly its Mott-Hubbard transition. In the future, we hope 
that it may be possible to use the multi-singlet insertion 
to resolve the question of the existence of a finite tem- 
perature Z2 breaking in this model and we are currently 
exploring this possibility. 

In summary, we have proposed and given evidence in 
support of the proposal that an operator which inserts a 
pair of singly occupied sites in a singlet spin configuration 
into a one dimensional Hubbard model acts as an order 
parameter for the Mott-Hubbard transition. The transi- 
tion occurs in this language because of the breaking of a 
spatial Z2 symmetry in the charge sector of the model, 
doubling the charge unit cell and causing an insulator. A 
natural conjecture based on this picture is that a similar 
Mott-Hubbard transition can occur in higher dimensions, 
involving the breaking of a charge sector Z2 symmetry, 
possibly at finite temperature. In support of this, we 
have given some numerical evidence for a broken Z2 for 
the two dimensional Heisenberg model and identified two 
experimental systems whose behavior is suggestive of the 
proposed, new kind of ordering. 
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